ASYMPTOTIC EXPANSION FOR THE WAVE FUNCTION IN A 
ONE-DIMENSIONAL MODEL OF INELASTIC INTERACTION 
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Abstract. We consider a two-body quantum system in dimension one composed by a test particle 
interacting with an harmonic oscihator placed at the position a > 0. At time zero the test particle 
is concentrated around the position Rq with average velocity ±vo while the oscillator is in its ground 
state. In a suitable scaling limit, corresponding for the test particle to a semi-classical regime with 
small energy exchange with the oscillator, we give a complete asymptotic expansion of the wave 
function of the system in both cases Rq < a and Rq > a. 



1. Introduction 



The analysis of the classical behavior emerging in a quantum system is a widely studied subject 
both from the theoretical and from the experimental point of view. It is generally accepted that such 
^ ' behavior cannot be explained by simply taking the limit ^ — )■ for the isolated system but the crucial 
Qs^ role of the environment must be taken into account. In particular the environment is responsible 
C^l for the suppression of the quantum interference between two different components of a superposition 
state of the system. The dynamical mechanism producing the suppression is known as decoherence 
f |BGJKS] . [GJKKSZj . for some rigorous resuhs see also [AFFlj . [UUF] . [D], [DuFTj . (DS]). 
In this context, a particularly important problem is the explanation of the appearance of straight 
tracks in a cloud chamber produced by an a-particle emitted in the form of a spherical wave. The 
problem was raised by Mott ([M]) who gave an heuristic explanation based on the analysis of the 
interaction of the a-particle with only two atoms of the vapour. In jPFTlj and |DFT2] the same 
problem was studied in dimension one and three, and a rigorous result was given up to second order 
in perturbation theory. In this paper we consider an even simpler situation, i.e. a test particle in 
dimension one interacting with one harmonic oscillator, and we want to give a complete asymptotic 
expansion of the wave function of the system valid at any order in perturbation theory. 
We want to emphasize that a complete and detailed analysis of this case is the starting point for the 
analysis of the more realistic model of an a-particle interacting with many atoms. 
The Schrodinger equation for the wave function \E'(t) of the system is 

= -^^Rm - ^^Mt) + ^rnu'ir - + XV{6-\R - r)Mt) (1-1) 
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In ( II .ip we have denoted by R, M the position coordinate and the mass of the test particle, and by 
r, m, oj the position coordinate, the mass and the frequency of the oscillator placed at the position 
a > 0. The smooth interaction potential is denoted by V and A > 0, 5 > are the coupling constant 
and the effective range of the interaction. 

We consider two possible initial state chosen in the product form 



. Mvn 



vi>±(i?, r) = — e±^^^r/(or-i(i? - R,)) ^,{r) (1.2) 

where cr, > Rq & rj & '5(M), \\rj\\L2 = 1, where S(R) is the Schwartz space. We have denoted 
by ±Vq, Rq the average velocity and position of the test particle at time zero and by ipo the ground 
state of the oscillator. More generally, eigenfunctions and eigenvalues of the oscillator are denoted 

by 



^nir) = ^Ml~\r~a)), nEN, 1 = \— (1-3) 

V rnoj 



En = nw[n + ^ (1.4) 
where 0„ is the Hermite function of order n. 

We want to study the asymptotic behavior of the wave function of the system under a suitably 
chosen scaling limit. More precisely we introduce a small parameter e > and consider a semi- 
classical scaling for the test particle for £ — )■ 0, i.e. 

h = e^ M = l a = e (1.5) 

For the physical parameters of the oscillator we fix 

m = e u = (1.6) 

which imply 7 = £, i.e. the oscillator is well localized around the position a, and En = e{n + l/2), i.e. 
the energy level spacing of the oscillator is much smaller than the kinetic energy of the test particle. 
Finally we choose the range of the interaction of the same order of the initial localization of test 
particle and oscillator and fix the coupling constant such that the deviation from the free dynamics 
is small, i.e. 

6 = e X = e^ (1.7) 
Under the above scaling, the system is described by the rescaled wave function ^'^(t) given by 

^f{t) = e-^^^^^^Jo (1-8) 
where is the Hamiltonian of the system 

H' = H'o+e^V' (1.9) 



1 

— + -(r - 
2 2^ ^ 



(1.10) 

V = V{e-\R-r)) (1.11) 
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and \1/^Q are the two possible initial states 

vi/Jo(i?,r)=V^±(/?)0^o(r) (1.12) 

i^fiR) = ^e^'^^''v{e-\R~Ro)) (1.13) 

cPKr) = ^Me-\r - a)) (1.14) 

The aim of this paper is to characterize the asymptotic expansion of the wave function of the system 
^^(t) for e^O. 

We list some notation which will be used in the following. 

a = {ai, . . . ,ai), da = dai . . . dai (1-15) 

I I I 

\a\ = '^^\aj\ ||a||^ = lojl^ a-b = '^^ajbj (1-16) 

j=i j=i j=i 

The norm of L^(M) will be denoted by || ■ \\lp and the inner product of L^(M) will be denoted by 
by (■, ■). The norm of the total Hilbert space L^(M^) will be simply denoted by || ■ ||. Moreover we 
introduce the norm 

k 

iii/iiu,p = Eii/'(-)'~'ii^'' (1-17) 

j=0 

where denotes the derivative of / of order j and (x) = (1 + x^Y^^. The symbol / indicates 
the Fourier transform of /. The propagator of the harmonic oscillator centered in the origin and 
corresponding toh = 'm = uj = l will be denoted by U{t). Its integral kernel is explicitly given by 

°° 12 2 

U{t; x,y) = Y: e-("+i)Vn(x)0„(y) = .-^^^t^^ (1.18) 

n=0 ^ ^ ^' 

Finally c will be a positive numerical constant whose value may change line by line. 
The paper is organized as follows: in section 2 we state our results and outline the strategy of the 
proof. The proofs are postponed to section 3 while in section 4 we explicitly compute the first terms 
of the asymptotic expansion and discuss a physical application. In the appendix we give the proof 
of proposition 12.11 



2. Results 

We first give some heuristic arguments in order to explain the qualitative behavior of the system. Let 
us consider the initial state "^^q, i.e. the oscillator in its ground state and the test particle initially 
described by a coherent state with average position Rq and positive velocity Vq. For small t > 
the effect of the interaction is negligible and then the state of the test particle will undergo a free 
semi-classical evolution. As the time increases, it is reasonable to expect two completely different 
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situations in the cases Rq < a (the test particle "hits" the oscillator) and Rq > a (the test particle 
does not "hit" the oscillator). 

In the case Rq < a, taking into account that the effective dimension of the oscillator and the range 
of the interaction are of order e, one has that the interaction will take place approximately at the 
(classical) impact time r given by 

r=\^^ (2.1) 

Due to the interaction, the state of the test particle and of the oscillator will be almost instantaneously 
modified. Then for t > r we shall find the oscillator in a superposition of stationary states while for 
the test particle we expect to find again a free semi-classical propagation, but of a modified coherent 
state. Therefore for this case the problem is reduced to compute the modification of the state of the 
system determined by the interaction at t — t. This will be done using a perturbative expansion of 
the wave function of the system and exploiting stationary phase methods to control each term of the 
expansion. Therefore we shall refer to this case as the stationary case. 

In the case Ro > a the situation is easier since we have only to verify that for any t > the state of 
the system has an essentially unperturbed evolution, with the oscillator in its ground state and the 
test particle evolving according to a free semi-classical propagation of its initial state. This will be 
done using standard non-stationary phase methods and therefore this case will be referred to as the 
non-stationary case. 

It is obviously true that if we consider the other possible initial state \E'^o then the opposite situation 
occurs, i.e. the case Rq < a corresponds to the non-stationary case while the case Rq > a corresponds 

to the stationary case. 

Let us outline the strategy of the proof for the characterization of the asymptotic behavior of "^^{t) 
(for \E'^ (t) the analysis is completely analogous) . 

The starting point is to represent the solution using Duhamel's formula 





After k + 1 iterations we get 



k+l 



^t{t) = e-^^^^ (^^+0 + E (^) j + (2-3) 
where we have denoted for any I — 1,2 . . . ,k + 1 

If{t) = i-iY fdsi . . . [dsi e'^^° V e-'^"° . . . eO"o e-'f^^°*+o (2-4) 
Jo Jo 

and the rest 7^|_|_i(t) is 

,{t) = -i fds e-'l^^V^e-'^^o Il^^{s) (2.5) 
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It is convenient to consider the expansion of the wave function "^^{t) on the basis of eigenfunctions 
of the oscillator, i.e. 

oo 

^tit;R,r) = Y,r^{t;R)^l{r) (2.6) 



n=0 



and 



From (12. 3p and (12.71) we have 



mR)= (0^,v[/+(t; /?,.)) (2.7) 

C30 

K{t)r = J2\\fnmh (2.8) 



n=0 



fc+1 



m = e-'^^'^ (5.0 e-^^^+ + e-^^("+^) ^^{<j>t,, , Jf (t))j + , 7^|(t)) (2.9) 

The next step is to give an explicit expression for the quantity (0^ , If(t)). We introduce the rescaled 
variable x given by 

X = e-\R- Ro) (2.10) 

and the following phase function 

(Pe{x} = —X - —Ro (2.11) 

e 

We have the following representation formula. 
Proposition 2.1. For any e > and n,l & N we have 

(0^,/f(t;£x + /?o,-)) = ^e*''^(^)(-2)' V 5n,ofdsJ\si^^...f\sJd$,Fr,{s,i-x)e^''-^'^^^ (2.12) 
where 

I 

Fr.{s,i-x) = r,{x-s-i)\[V{Q<K^n,{iM^'-^^^^^^^^ (2.13) 

i=i 

I I 

^r^{s, $.)=Y^ $,(Sj, 0) = 5^ [(« - ^0)0 + (%+l - - ^OSjO] (2.14) 

i=i i=i 

anc? n^+i = ,^0 = 0. 



The proof is a long but straightforward computation and it is postponed to the appendix. From 
formula (I2.12p one sees that the problem is reduced to the analysis of a highly oscillatory integral. 
The asymptotic behavior for e — )■ is then characterized by the presence (or the absence) of critical 
points of the phase in the integration region. The unique critical point of the phase is easily 
found 
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s'' = sgn{a- Ro)t, t = (r, . . . , r) (2.15) 



r = (ej,...,en, q = ^^^^^^, j = i,...,/ (2.16) 

For Rq > a the critical point lies outside the integration region and then we can use standard non- 
stationary phase arguments to show that the contribution of the integral can be made arbitrarily 
small. The result is summarized in the following theorem. 

Theorem 2.2. Let us fix k t> and assume that 1 1 |V^| | < oo. 
If Rq > a then 



and if Rq < a then 

The remainders M^{t) satisfy 

where Dk depends on t, vq, V , rj. 



^tit) = e-'^"^^t,o+KHt) (2.17) 
vl/;(t) = e-'^^'^^^.+Kit) (2.18) 
\\Aftm<D,e'^' (2.19) 



For Rq < a the critical point lies on a one dimensional subset of the boundary of the integration 
region. This means that the explicit computation of the asymptotic expansion in power of e is rather 
involved (see e.g. [BHj.[F].jH]). In order to simplify the analysis, we shall exploit the fact that the 
phase in formula fl2.12p is linear in the variable More precisely, we have 

y"d|F„(s,^;x)e-^*"("'^) = e^''^^^-' j d^Fr,{s, x)e>^'-^>^ (2.20) 

where we have used the notation introduced in fl2.15p . fl2.16p . Exploiting fl2.20p in fl2.12p . introducing 
the new integration variables Zj = e'^vq^sj — t), j = 1, . . . , / , and defining 

fie = {;2 G I - e'^voT < zi < Z2 < . . . < zi < e'^voit - r)} (2.21) 

we get 

(0^,Jf(t;£a; + i?o,-))=f— ) -^e'^^^^^ e^''^ T 6n,o [ dze'^''-^ [d^Fr,{T + ev,'z,C,x)e-'^-^ (2.22) 

Equation 02.221) suggests that in order to have an asymptotic expansion of If{t) it is sufficient to 
consider the Taylor expansion of Fn{T + ) and to prove that for each term we can extend 

the integration w.r.t. z to 

fio = e M' I - oo < zi < Z2 < . . . < zi < oo} (2.23) 



paying only a small error. Let us rewrite {(j)^^, If {t)) in a more convenient form. From fl2.22p and 
(12.131) we have 



mt;ex + Ro,-))- 



^ 0/ V ni...n;=0 -^^"^ j=l 

We introduce the following matrix valued function 



M{z) 



Z2 Z2 2:3 
^3 ^3 2:3 



(2.25) 



\zi Zl Zl ... Z;/ 



and we write 



Furthermore we have 



i=i V fc=i / 



(2.26) 



(2.27) 



The series over tt-i . . . n; = can be explicitly computed and each of them reconstruct the propagator 
of the harmonic oscillator. In fact 



oo I 



ni...n;=0 jr=l 



riiO 



ni...?i;=0 3 = \ 



8n 



lO 



= (2vr)-^-/2,-(n+l/2)., (0„,C.(O)n^(O) 

where in f l2.28p we have introduced the notation 

G(0=e-'«"-[/ 



(2.28) 



^0 



(2.29) 
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In the end we arrive at 



, If (t; ex + i?o, ■)) = f-^^l ^ e^^^(^)e^"^ / dz e^"^' / e-^-^e^K^^^i «^)'(0„, C,(|)) 



n ^(0) - r ■ I - ^i^o'^z ■ ^)e^^-'^'(")^ (2.30) 



Our main result on the asymptotic behavior of for e — in the stationary case is summarized 
in the following theorem. 

Theorem 2.3. Let us fix k eN, t > r and assume that \ \\V (■)^| | |a:+2,i < oo. 
If Rq < a then 

^tit) = e-^l^^o (^-^e^'^^g.^X+j +5+(t) (2.31) 
%{t) = e-^^^o (^-^e^^^Y.e'^I,^ +S;{t) (2.32) 



and if Rq > a then 



where the coefficient of the expansions are given by 

h 

I±=#o X,± = 5^Xf'^ h^l (2.33) 

1=1 

and 



■'ni 



I h—l / . \ h—l—m 

n ^(^.O ^"^'^ E ^'"(^ - ^ ■ ^) (-^ ■ 0™ ( ■ M{z)^ ) (2.34) 



j = l ni=0 

The remainders iS^ (t) satisfy 

\\S^{t)\\<Cks'^' (2.35) 

where depends on t, vq, V , rj. 



Theorem 12.31 shows that for t > r we have again a free evolution of a modified initial state. 

We notice that the factor e"^/^ in f l2.3ip . fl2.32p takes into account the change of scale between R 

and x and therefore £~^/^||X;'''^|| does not depend on e. 



3. Proofs 



In this section we shall prove theorems 12.21 and 12.31 For the sake of simplicity we shall drop e as index. 
We are not interested in explicitly determining the dependence on k of Ck and and therefore many 
combinatorial factors appearing in the following formulas will be included in c. Let us start with the 
non stationary case. 

Proof of theorem 2.2 

We shall give the proof only for '^^{t) with Rq > a, since the case of with Rq < a can be 

treated exactly in the same way. Let us consider the expansion (12.31) with k = 



Taking into account (12.51) . we have 

||^+(t)_e-f ^6^+11 ^ ^i + t\\V\\Lo.) sup (3.2) 

From (13. 2p we see that the problem is reduced to estimate ||/i(t)||. Formula (I2.12p for / = 1 reads 



(3.1) 



hit; ex + i?o, ■)) = ^ e^^^(^)ei"-y^tis j ViO^oiOvi^ - sOe^^«'e-"«e-'^H-«-^^+^«; 

(3.3) 

Therefore we have 



n=0 
oo 



ds di ^(O0n0o(O^(^ - sOe-^'^ e-'-^e 



dx / dsdsld^dC ViOi^o{OVi^')i^oi^')vix- sOvix- s'Oe'^'^'e-'r^'^" 

n=0 J Jo J 



The series over n in (13. 4p can be explicitly computed as follows. 

oo ^ -l oo « 

E^(0£S(e)e-^("+^/')'^ = 7rJ2 dydy'Uy)My)e-''^Uy')Myy'^'e 



n=0 



n=0 



/w .. n(s'-s) 



Zn \ 



s — s 



(3.5) 
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Substituting f l3.5p into fl3.4p we obtain 

The estimate of fl3.6p can be easily obtained exploiting a non-stationary phase argument. We consider 
the identity 

^-IMr+s)^ = dk+le-iMr+s)^ (3.7) 

and we integrate by parts k + 1 times. Hence 

.2k+2 p pt 1 r ^=+1 /; , i\ ^+1 



bl='=+i |j'|=fc+i 



■ r]^'^{x - s'O {-sy'^V^'^{i')&>le^^''^'\-'''^' {-ix'y'^ (3.8) 

where 

(3,9) 



Jl! J2! Js! J4! J5! 

is the standard multinomial factor. From fl3.8p we have 



lil=fc+l \3'\=k+l 

■ \v'Km^\'M''\^i&^^WsAi.O\WKx - sa\ s''^\vHO\\si'\^x'Y'^ (3.10) 

We interchange the integrals and we estimate the integration w.r.t. the variable x using Cauchy- 
Schwartz inequality 



^ dxW'{x - sOM^W^{x - s'OWx'y^ ^c\\r]''{-y'\W{s^Y'\\r]'^{-y^W{siy^ (3.11) 
Moreover from the definition (13. 5p and Cauchy-Schwartz inequality we have 



^ \W'<Po\\\W'^<Po 
1 



^ rfj5 + ^)r(ju^) (3.12) 
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If we substitute ( 13. lip and (13.121) into (13.101) we have 

fc+i i>t f '^+1 



iii=fc+i 



/f\k+2 ^=+1 



k+l 



^ ce 



.+1 {tr' 



jl,...,j5=0 
|i|=fc+l 



1^1 1 |fc+l,2| 1 1^1 1 |fc+l,l 



The proof of the theorem now follows from (13.21) and (13.131) . 

For the proof of theorem 2.3 it will be useful the following technical lemma. 



(3.13) 



□ 



Lemma 3.1. Let Cti...t„_i : IR" ^ L^(^) be defined by 

C*i...*„-.(0 = e-«-t/(t._i)e-«"- 



^7(ti)e-^«i-0o 



(3.14) 



Then Cti...t„_i ^ C°°(M"; L^(R)) and for every multi-index ct = («!,...,«„) G there exists Ca, 
independent ofti for i = 1, . . . ,n — 1, such that the following estimate holds 



i2 \h=l 



(3.15) 



Proof 

The derivative d^(ti...t„-i{^) is explicitly given by 

dtCt^...t^-^i^) = (-2)l"le-^f"-Q""[/(t„_i)e-^«"-i-Q°"-i---f/(ti)e-^«-g"i0o 



(3.16) 



where Q denotes the multiplication operator by the independent variable y. In order to estimate 
the norm of (13. 16^ we shall move all the Q's to the right until they act on 0o using the following 
identities 



Q U{t) = U{t) (Q cos t - D sin t) 
D U{t) = U{t) {Q sint + D cos t) 



(3.17) 
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where D = i-^. Let us first see how the strategy works for the case \cx\ = 1. 



-i)e"*^"' ■ ■ ■ U{tj^i)e-'^'-^- (Qcostj_i - L)sintj_i - sint^^i) U{tj_2) ■ ■ ■ U{ti)e 



Iterating the procedure we arrive at 



(3.18) 



di^jCti...t„-i{i) 



-tie 



i-1 i-1 i-i i-i 

' f/(ti)e-'«i- COS (J2 ^ - ^ sin (Y1 ^ " 5Z ^™ sin t, 



m=l 



m=l 



m=l 



h=m 



h (Po 



(3.19) 



From equation (13.191) we obtain the estimate 

m=l \m=l 

Derivatives of any order can be estimated by iteration of the procedure used in (I3.19p . starting from 
the most left derivative and going to the right. This leads to 

/i-i \ /i-i \ j-i /j-i 



n 



m=l 



(3.20) 



L2 



From estimate (I3.20p it is clear that there exists Cq, depending on suitable norms of 00, such that 
(Km holds. 



□ 



In the next proposition we derive an expansion for the quantity Ii{t)) which is the key ingredient 
for the proof of theorem 2.3. 

Proposition 3.2. Let us consider ^'^q as initial state and let us fix Rq < a, t > t and /c, / G N, with 
I ^ k. Then there exists a constant C, depending on k,l,vo,T,t,r],V , such that 

k 

{<Pn,m) = 5^e^0„,Xf'+) +e'=+^(0„,r(t)) (3.21) 



h=l 



where 



,Xf-+(ex+i?o,-)) = ^ ^e-^-(-)ei-/c^zet("+V2). r d^,-^t^e^r(J:U,,Y ^^^^^^^^^^ 



I h—l y . \h—l—m 

n ^(^.) ^"^^^'^ E ^"(^ - ^■^) (-^■^)'" ( \^ ■ ^^(^)^ ) (3.22) 



m=0 



and 
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(3.23) 



Proof 

As a first step we consider the Taylor expansion of tlie e-dependent part of the integrand in (12.301) 
up to order k — I. We have 



j=0 m=0 ^ ^ 



m=0 



.ei^«-n-i?(l_^)'c-*+i (324) 
The corresponding expansion of {(pnyliit)) is given by 

k 

{<Pn,m) = Y,^'{<PnJiW+e''-'{<Pn,m) (3.25) 
h=l 

with 

/i!(zV27r)'fg^ Jf7, J 

■ n ^(0) e"'^^" E ^"(^ - ^■^) (-^■^)"^ l^-m^)^ (3.26) 



and 

k-l+l 



m=0 



(A;-/ + l)!.r'+^..^ 

rf0r7™(x - T-C - £ei;o"'^-0 (-^-0" U ^•^(^)^J e^^-'^'^"^^(l - 0)^^-'+^ (3.27) 



We have to prove that each term in the r.h.s. of (I3.25P is well defined, i.e. (t) and T{t) must be 
elements of L^(M^). We show this fact only for the most singular term T'(t), since the other terms 
can be treated in a similar way. 
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The square of the norm of T^{t) is given by 



/ . \ k-l+l-m . ^1 

(\ fe— i+l— m , 
^' ■ Miz')n e-^^'-^'^^'^^a - e')'-'+' (3.28) 



where, once again, we have reconstructed the propagator of the harmonic oscillator exploiting the 
formula 

CA^'),u(^^^ = get("+V2)(.,-4)^^^^^^(^)^^^^^^^,(^.)^ (3.29) 

In order to control the convergence of the integrals over fl^, we integrate by parts k + 2 times w.r.t. 
to and to respectively, for j = 1, . . . ,1. For the sake of simplicity, we explicitly perform the 
computation in the case I — k — 2, the general case being more delicate only from the notational 
point of view. Then we are reduced to estimate the following quantity for m — 0,1 



■ ViCi)ViC2)H'i)V{Q e-(^i+^^)-e^(«i+^^)- d0r)^{x - r • £ - e0v^'z ■ ^) {-z ■ 

Jo 

• U ^ • M(z)d e^^-^(^)^(l - 9) de'rj^ix - r • - ee\'z' ■ {-z' ■ 

= / dzdz'G{z,z') (3.30) 

We shall derive a pointwise estimate of G. Let us introduce the shorthand notation dl" = and 
di' — dp and let us consider the identity 



-tZi 



)4g-i|.^ = d^e-'^"^ (3.31) 
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Then integrating by parts four times w.r.t. and ^[ in f l3.30p we have 
zt4G{z,z')= [dx[ddd^'e-^^-e^^'-'V{^2)V fl (^) E 

ii,...,i7=0 ^-^^ j[,...,j',=0 



^2 ~ ^2 

^^0 



b1=4 |j'|=4 



1— m 



(1 - e) di'r]"'{x - T-^ - eOvq^z-^,) di' I (-Z ■ $,r [2 ^ ' ^^^^^ 



1— m 



(3.32) 



Using fl3.15p and exploiting the fact that \\evQ ^z\\ ^ \t\ for z G fie, we have 

z^2'i^|G'(z,z')| ^c|2||;2'| jdx j^dode' {i-e){i-e') jd^d^'\v{^2)\\v{Q\ 



(3.33) 



h,-,j7=0j[,...,j!^=0 
I J 1=4 1/1=4 



We use Cauchy-Schwartz inequahty to estimate the integral w.r.t. x and we obtain 



ztz[^\Giz,z')\ ^ c(t)*|z||2" 



^ c{tf\z\\z'\ 



il J2 J3J4=0 

^ii+i2+i3+i4=4 



jlj2=0 

.ii+i2=4 



ii,j2=,o 
ii+i2=4 



^c|;2||;2'|(t)«|||r^-|||l,|||\/(-)1||t, 



(3.34) 



Since the variables zi, z'l play no special role, we can repeat the same computation for each variable 
and we still obtain the same estimate as fl3.34p . Summing up we arrive at 



|G(^,z')|^c(z)-^(z')-^(t)«|||r/-|||l2lim)1llli 



(3.35) 
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The above estimate guarantees that (see f l3.30l) 

E Ell'^-WlliO ^c(^)'lll^(-)'lll4,i E lll^'^llk^ (3.36) 

r!i=0,l \n=0 / m=0,l 

in the particular case I = k = 2. In the general case, exploiting the same argument, one obtain 

\G{z, z') I ^ c {t)'^'^'^ 1 1 |r^-| 1 1 ^2,2 I \\V{-r\ 1 1^^2,1 (3.37) 

and therefore 

k-l+l / oo \ 1/2 k-l+1 

wmw^ E Eii'CWiiiO ^c(^)'^'iii^(-)'iiiUiE 1 1 1^'"! 1 1^^+2,2 (3.38) 

m=0 \n=0 / m=0 

Using the same kind of arguments, it is easily shown that also (t) belongs to L^(]R^) and this means 
that equation (13.251) is well defined. To conclude the proof of the proposition it remains to estimate 
the difference Xf'+ - (t). We have 



I h—l / . \ h—l—m 

■ {^n, cm ■ n ^(^.) ^"^'^ E ^'"(^ -r-^){-z- \i ■ M{z)i (3.39) 



(3.40) 



j=l m=0 

If we repeat once again the integration by parts procedure outlined above we obtain 

||Xf'+-/f(t)Kc(t)'=+2|||\/(-)^|||Ui E lll^'^IIU+2,2 / dz 
We observe that 

Therefore the expansion (13.2ip is proved with 

k 

m = fit) + E^'-'"'(^'^ - Iiit)) (3.42) 

h=l 



□ 



We are now in position to prove theorem 2.3. 
Proof of theorem 2.3 
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We give the detail for asymptotic expansion associated to \E'^q with _Ro < cl- The case with \I/^q with 
Rq > ais completely analogous. We consider formula fl2.3p and we expand each Ii{t), for / = 1, . . . , fc, 
according to proposition 13.21 Thus we obtain f l2.31|) f l2.33|) and f l2.34p . It remains to prove that the 
rest is of order e'^'^^. From (12. 5p we have 

ll^fe+iWII ^ II^IU- / ds \\h+i{s)\\ (3.43) 
Jo 

It is straightforward to notice that, under the assumptions of theorem 12.21 one has the estimate 
||/fc+i(s)|| ^ ce^~^^ (we just have to proceed as in proposition 13.21 integrating by parts k + 2 times 
without the Taylor expansion) and therefore the proof is concluded. 



□ 



4. An application 



Here we give an explicit computation of the first and second order terms of the asymptotic expansion 
of ^E'^(t) in the stationary case and discuss an application. From theorem 2.3 we have 



^ e'^^ (Xo+ + ell'^ + (X^+ + 1',^)) 



+ 0(r 



(4.1) 



The first order term is given by 

1 



,X^+(£X + i?o,-)) 



1+1/2 . 



iVoV^TT 

The integral over z can be explicitly computed 



(4.2) 



i(!i±iZH_5), /n + 1/2 
dze^ "0 ^' = 2'no | — 



and then one obtains 
(0„,X^+(£x + i?o,-)) -- 



O-TT ^ » (n+1/2)^ 

ee e 



; n + l/2 



6„,e "0 '(t)o)V 



n + 1/2 



J 1 — L — '; 



7] \X — T- 



(4.3) 
_n+ 1/2 



^^0 



(4.4) 

We also compute the the two terms appearing at the second order. This calculation shows a glimpse 
of the difficulties appearing in the general term of the asymptotic expansion which are not hinted by 
the first order term. The computation is a bit more delicate and it involves 5, 6' and principal value 
distributions. The first one is the following 



,Xf'+(£X + i?o,-)) 



i 



(4.5) 
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Also in this case the integral over z can be computed explicitly. 



/ 



A K- — ■ 



^^0 



Therefore 



1, 



n + 



-e= < IT- 



1 



2 \ — * X 



n 



■q \ x — T 



Vq J 2 vl 
e ' -0 ■ . (/)o)\/( 



^^0 

e "0 



?7 I a; — T- 



^^0 



^0 



I X — T 



Tj \X — T 



''0 



T] [X — T- 



n + 



vo J 2 



2 / , 

T] \ X — T 



n + 



2^ ^3 



.n+A 

I ^: 

"0 



77 I X — T 



Vo J 2 



Vo 



7] [ X — T- 



Vo 



,n + 



+ e^' "0 {c|)r^,e-'^■c|>o)V{^-^^){-ix)e 

Vo 



. n+l 
-t -X 



n + h\ n + „ 
T] \ X — T I t 



(4.6) 



^^0 



Vo 



- r 



2 vl 



-\- l \ T]' [ X — T 



^^0 



+ T- 



n + i 



2 



^^0 



T] \X — T 



1^0 



(4.7) 



The second term contains a less trivial integral over z. 



+1/2 . 



(0„,X|'+(£x + i?o,-)) = ^TT^e""" / dz,dz2e'^'' d^ic^^ e-*(^i«i+^^«^)e^^«i+«^) 

7Z7;o27r 



e-^«^-t/ ( ^ ) e-^«-0o)t^(6)t^(e2)e^^(«^+«^^^r;(a: - T(ei + 6)) (4.8) 



If we expand the harmonic oscillator propagator 



(4.9) 
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we are reduced to compute 



(4.10) 

We have 

"°°d.,e-*"(^'«-^)= 2rf f6 + 6 + (4.11) 



OO 



^^0 



We already remarked that the stationary phase expansion is different with respect to the usual one 
due to since the stationary point lies on a one dimensional subset of the integration domain boundary. 
The appearance of distributions, like the Principal Value in fl4.12p . is a manifestation of this fact. 
We have 



_ _ !!i±i±Z! V 1 pv 



vo J 2m ^2 - ^1 



rra+l+n 



We notice that, following the same kind of arguments, in principle one can explicitly compute also 
any other higher order term of the asymptotic expansion of '^^{t). 

We want to conclude with a brief discussion of a possible physical application of our result. Let us 
consider the following initial state 

vi>,,o(i?,r) =a, I -^r^{e~\R-R^)) + —^r^{e-\R + R^)) ] ct>l{r) (4.13) 



where f o, -Ro > 0, -Rq < a and is the normalization factor 

-1/2 



ae = ^(^ + Jdxr]{x)r]{x -2Roe~'^)cos2ivoe''^{x - Roe ^)^ (4.14) 

We notice that by repeated integration by parts one easily obtains = ■^ + 0{e''), for any k eN. In 
f l4.13p we are assuming that at time zero the test particle is described by a coherent superposition of 
two wave packets with opposite average momenta and the oscillator is in its ground state. Exploiting 
theorems 2.2, 2.3, we explicitly write the first order approximation of the time evolution of (14.131) 
for a fixed t > r 

vi/,(t) = e-^^^oe, + ^,(t) (4.15) 
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where < ce^ and 



e,(i?, r) = a, e-2^^^« e„,,(i?) </)^(r) (4.16) 

n=0 

i>g _ n + 1/2 ^ ^ / _ 1 



+ /3n,.v^e^^ ^rry(£-i(i?-i?o)-r(n + l/2)t;o-i) (4.17) 

\/2tT i /" (?i+l/2)iio , \ j ("+l/2)^ .n+1/2 ^ 

/9n,. = ^eH^^+^^Je ^^(0„,e-'^-0o)^K + 1/2)) (4.18) 
It is also convenient to consider the Fourier transform Qn,£ of 'S)n,ei which is exphcitly given by 

+ /3„,,£3/2g-<A-^+^)(^^+iJo)^^^^_^-i^^^^-i(^^^/2)) (4.19) 

Starting from fl4.15p . our aim is to compute approximate expressions for the probabihties of the 
outcomes of the measurement of some interesting observables relative to the system which could 
possibly be performed for t > t . In particular we are interested in: 

- the probability to find the energy Eq or Ei for the oscillator, denoted by Vq, Vi respectively; 

- the probability to find a positive momentum for the test particle and the energy Eq or Ei for the 
oscillator, denoted by P+_o and V+^i respectively. 

In order to control the approximations, we observe that for any projection operator P in L^(M^), 
which commutes with H^, one has 

|(vl/,(t),PM/,(t)) - (e^PGe)! < ||f,(t)||(2||P0,|| + \\SM\) < ce^(\\Pee\\+e^^ (4.20) 

From (14.151) . a direct application of the Born's rule and (14.201) we obtain 

Vo = \a,\^ jdR \QoAR)\' + Oi^') = 1 + Oi^) (4-21) 
Pi = \a,\' jdR\Q,,,{R)\' + 0{e') = ^\f3,J' + 0{e') (4.22) 

Moreover 



|2 / jr^lS /r^M2 , /n^^2\ ^ 



V+^o = WeV j dK\Q,^,{K)\' + 0{e') = -+0{e) (4.23) 
= dK \ei,eiK)\' + Oie') = - + Oie') (4.24) 
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Hence we have 

^ +Oie) (4.25) 
1 + 0{e) (4.26) 



Vo 2 



Pi 

The Lh.s. of formulas fl4.25p . fl4.26p could be interpreted as a sort of conditional probabilities. More 
precisely, (I4.25P says that if the oscillator remains in the ground state then the probability to find 
a positive momentum for the test particle is approximately one-half. This essentially means that in 
this case the superposition state for the test particle survives. On the other hand, (14.260 says that 
if the oscillator is in the first excited state then the probability to find a positive momentum for the 
test particle is approximately one. In such case the test particle is described by a wave packet with 
positive momentum moving to the right of the oscillator. It should be stressed that both "histories" 
are contained in the complete wave function of the system and a choice is made only when a real 
measurement of the energy of the oscillator is performed. 



5. Appendix 

In this appendix we give the proof of proposition 12.11 which is essentially based on algebraic manip- 
ulations and the use of Fourier transform. From (12.41) we have 



„.-n JO JO Jo 



ni...n;=0 



e^^^^T/e-^'^^V^,.,) . . . (0„„e^ii^«\/e-^ii^V„,)^+ (5.1) 



For the generic scalar product appearing in the integrand in (15. ip we have 

e'^^^V^e-'^^'V^) = e^("-'")^e'^'^°V;me"'^''" (5.2) 



■"nj ' 



where 



V:miR) = l I {e^\r-a))V{e-\R-r))<Pm {e-\r - a)) 
^dx 0„(x)0m(a^)V" [e~^{R — a) — x) 

^d^ViO^miOe-'^"^ (5.3) 
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Furthermore for any g G we have 



)= Jdke''^e''^'e-''^i'-'lTg{k + 6-'^) 
e-''^ fdke'^^e-'''^^g{k + e-^i) 



\/27r 

= e''^^ e-'^^g{R-esO (5.4) 
From (|521), O, (El]) we conclude 

Substituting (15. 5p into (15. ip and taking into account the exphcit expression of ip'^ we obtain (I2.12p . 

□ 
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